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Periodic waves are investigated in a system composed of a Kuramoto - Sivasliinsky 
- Korteweg - de Vries (KS-KdV) equation linearly coupled to an extra linear dissipative 
one. The model describes, e.g., a two- layer liquid film flowing down an inclined plane. It 
has been recently shown that the system supports stable solitary pulses. We demonstrate 
that a perturbation analysis, based on the balance equation for the net field momentum, 
predicts the existence of stable cnoidal waves (CnWs) in the same system. It is found 
that the mean value uq of the wave field u in the main subsystem, but not the mean value 
of the extra field, affects the stability of the periodic waves. Three different areas can 
be distinguished inside the stability region in the parameter plane (i^,Mo), where L is the 
wave's period. In these areas, stable are, respectively, CnWs with positive velocity, con- 
stant solutions, and CnWs with negative velocity. Multistability, i.e., the coexistence of 
several attractors, including the waves with several maxima per period, appears at large 
value of L. The analytical predictions are completely confirmed by direct simulations. 
Stable waves are also found numerically in the limit of vanishing dispersion, when the 
KS-KdV equation goes over into the KS one. 
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§1. Introduction 

It is commonly known tliat periodic waves, as well as their solitary-wave (SW) counter- 
parts, are coherent structures of fundamental importance in the study of physical systems 
governed by nonlinear evolution equations Stability of these solutions under small 
perturbations is an issue of basic interest in its own right and for applications, that is why 
it has been investigated in many works (see, e.g., refs. 0-0 and references therein). In 
particular, the stability of periodic waves in perturbed Korteweg - de Vries (KdV) and 



Benjamin-Ono equations was investigated, respectively, in refs. p-[|Tl|] and [0, [U 

An important one- dimensional wave-generating model that combines conservative and 
dissipative effects is a mixed Kuramoto - Sivashinsky (KS) - KdV equation, 

'^t ~l~ UUx ~\~ Uxxx CH'^xx I'^xxxx i (1) 



with a;,7 > 0, which was first introduced by Benney |]T^ and is often called the Benney 
equation. This equation finds various applications in plasma physics, hydrodynamics and 
other fields |TH]-[|T^. Localized structures are, obviously, important objects in systems 



of this type lHj-l^; however, SWs cannot be stable in the Benney equation proper, as 



the zero solution, which is a background on top of which SWs are to be found, is linearly 
unstable in this equation due to the presence of the linear gain, which is accounted for by 
the coefficient a in eq. (|l]) 0. A stabilized version of the Benney equation was recently 
proposed by the present authors in ref. ||22|. It is based on the KS-KdV equation for a 



wave field t), which is linearly coupled to an additional linear dissipative equation for 
an extra field f(x, t), that provides for the stabilization of the zero background: 

~\~ UUx ~\~ Uxxx Vx OlUxx 'yUxxxxy (2) 

Vt + CVx - Ux = TVxx- (3) 

Here, a, 7 and F are positive coefficients accounting for the gain and loss in the u- 
subsystem and loss in the v - subsystem, respectively, and c is a group-velocity mismatch 
between the fields. The system describes, for instance, the propagation of surface waves in 
a two-layer liquid film in the case when one layer is dominated by viscosity (see a review 
^ and ref. |^ and references therein). 

In the work it was shown that the stability of the zero solution to eqs. (0) 
and (H) is compatible with the existence of stable SWs. Actually, SW solutions for this 
system were obtained, with a high accuracy, in an approximate analytical form by means 
of a perturbation theory based on a balance equation for the net wave momentum. The 
existence and stability of these solitary pulses was confirmed by direct simulations. In 
fact, it is the first example of stable SWs in models of the KS type. 

In contrast to SWs, periodic waves can be stable in the Benney equation [p|-P|,[p^ 



25 1 . Periodic waves are more likely than SWs to be stable in dissipative nonlinear systems 
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because the stability of the zero solution is not a necessary condition in this case. Further- 
more, in experiments it is easier to observe periodic waves than SWs, this is why detailed 
study of their existence and stability is a relevant problem. 

A purpose of this work is to find periodic-wave solutions to the system of eqs. @ 
and and to investigate their stability. In §2, we find a family of exact cnoidal-wave 
(CnW) solutions to a zero-order system, in which the gain and loss are dropped. In §3, 
a perturbation theory for the periodic waves is developed by treating the gain and loss 
as small perturbations. To this end, the balance equation for the net field momentum is 



employed, following the technique developed for SWs in ref . . We demonstrate that a 
condition for equilibrium between the gain and loss may select two (or none) particular 
steady-state CnW solutions out of the continuous family existing in the zero-order system. 
In the case when two solutions are selected, the one with a larger amplitude is stable, 
while the other one is unstable, playing the role of a separatrix bordering an attraction 
domain of the stable CnW. 

In §4, we present results of direct numerical simulations of the full system, which turn 
out to be in very good agreement with the analytical predictions. In §5, we briefly con- 
sider the special case of the system with the vanishing dispersion, i.e., with the KS-KdV 
equation (|^) replaced by the KS one. In this case, direct simulations demonstrate the 
existence of stable periodic waves. The paper is concluded by §6. 

§2. Cnoidal Waves in the Zero-Order System 

In this work, we focus on the study of steady traveling-wave solutions to eqs. (0) and 
(§) in a spatially periodic domain < x < L with boundary conditions (b.c.) 

u{0) = u{L), u^{0) = u^{L), M^^(O) = u^xx{0) = u^^^{L)] 

viO) = v{L), vM = Vx{I^)- (4) 

Due to the imposed periodic b.c, the stability results obtained below (in §4) by means 
of direct simulations are limited to perturbations which have the same periodicity as 
the unperturbed solution. In a long system, a periodic wave may be subject to extra 
instabilities against periodicity-breaking perturbations. However, in this work we will 
display generic examples of stable waves with a fundamental period L/N rather than L, 
with N taking values, typically, 3 or 4. Those examples will demonstrate that, as a matter 
of fact, the waves may be stable against periodicity-breaking perturbations too. 

The group- velocity parameter c may be set equal to zero in eq. (||), as a substitution 

U ^ U + Uq, X ^ X + UqI (5) 

transforms eqs. (|^) and (^ into the same system, with c replaced by c — mq. In other 
words, we can eliminate c by subtracting the constant (background) value uq = c from 
the field u{x,t) (obviously, this transformation could not be applied in the case of SWs, 
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for which the zero value of the background must be kept). Therefore, c is assumed to be 
zero hereafter. 

In the zero-order approximation, corresponding to a = 7 = F = 0, eqs. (Q) and @ 
take the form of a conservative system consisting of the KdV equation coupled to an extra 
linear one: 

Ut + uu^ + u^^^ = v^, Vt = (6) 

(recall we have set c = 0). In ref. it was demonstrated that the conservative system 
(P) is not an integrable one. 

It is well known that the KdV equation proper, 

Ut + UUx + Uxxx = 0, (7) 

admits a family of periodic-wave solutions that have an arbitrary period L and travel at 
an arbitrary velocity s. These solutions are often called cnoidal waves (CnWs), which 
stems from the symbol cn for the elliptic cosine. Making use of definitions 

^ = x-st, r]L = 2K{m), (8) 

CnW solutions to the KdV equation can be cast in the form 

u{x,t) = I27f [dn2(r/^,m) - E{m)/K{m)\ , (9) 

s = Arf [(2 - m) - 3E(m) /K{m)] . (10) 

Here, t] is an intrinsic parameter of the solution family, which determines the amplitude 
and shape of the wave, dn(^, m) is the Jacobi's elliptic function with the modulus q = i/m? 
and K{m) and E{m) are complete elliptic integrals of the first and second kinds. In the 
limit of m — i> 1, the CnW solution goes over into the KdV soliton. Note that the solution 
(D has a zero mean value of m(0- 

On the basis of the CnW solutions to the KdV equation presented above, it is easy to 
generate a family of periodic- wave solutions to the zero-order system (^, making use of 
the same definitions (||) as above. The result is 

u{x, t) =Uo + I27f [dn^(r7^, m) - E{ra) /K{m) \ , (11) 
2 

vix.t) =Vo- 12— \dn^(ri^,m) - E(m)/K(m)] , (12) 

S L -I 

where the wave's velocity s is related to the mean value Uq of the u - field as follows: 

uo + V [(2 - m) - 3E{m)/K{m)] =s-^. (13) 
Lastly, f is the mean value of the w-field. 
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^3. The Perturbation Theory for Cnoidal Waves 



In this section, we perform a perturbation analysis of the full system of eqs. (H) and 
(I) in the presence of the gain and loss terms. First, we notice that eqs. @ and (||) 
conserve two "masses", 

M = u{x) dx, N = v{x) dx 
Jo Jo 

(which indeed have the meaning of masses in the application to liquid-film flows, see refs. 
0, and references therein). The net fleld momentum P for periodic solutions is deflned 

as 



P 



L r 



u^{x) + v'^(x) 



dx. 



(14) 



Unlike the masses, the momentum is not conserved by eqs. (0) and (§), but it is conserved 
in the zero-order approximation, setting a = 7 = F = 0, therefore the balance equation 
for P may be naturally used as a basis for the simplest perturbation theory 

With regard to the periodic boundary conditions (^, an exact evolution (balance) 
equation for P can be easily derived: 

^ = ^ {aul - jul^ - Tvl}j dx. (15) 

In the first-order approximation of the perturbation theory, we assume that the system 
inherits the unperturbed CnW solutions in the form of eqs. (pH]) and ([T^), while the 
amplitude parameter t] may vary as a function of the slow time, T = et. In particular, 
the momentum corresponding to the unperturbed solution is 



Po 



1/ 2 



1 + 4 



(16) 



2 ^ \ 

Calculating the integral on the r.h.s. (right-hand side) of eq. (plSj) analytically and sub- 
stituting the expression (O) in l.h.s., we derive an evolution equation for the amplitude, 



drj 



6r]^s^ {s^ + 1) 



105co 
where we define 

Co = 



3(s2 + l) 



4(s2 



8070377^ + 28ci (a - r/s' 



(17) 



(m - l)(3m - 2)ir(m) - 2(1 - 2m)E(m) 
-3 [E{m) + (m - l)K{m)f /K{m), 
C\ = {m — l){m — 2)K{m) — 2{1 — m + m'^)E{m), 

C2 = {m-l){m'^ + 2m-2)K{'m) -{2-3m-3m'^ + 2m^)E{m). 

In the subsequent analysis, it is more convenient to follow the evolution of the wave's 
velocity s, rather than the amplitude 1]. To this end, we express rj in terms if s, making 
use of eq. (|I3D: 



S — 1/s — Uq 



(2-m) -3E{m)/K{m) 



(19) 
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a corollary of which is 



ds _ 8t] [(2 - m) - 3E{m)/K{m)] 



Using eqs. (|T9|) and (0), we derive, after some algebra, the following evolution equation 
for s from eq. ([T7|): 



ds Cs^{s — 1/s — uq) 



where 



dT 3 (s2 + 1)' - 4 (s2 - uos - 1) 
47C2 



(s^ + aas^-s-ao), (21) 



7co[(2 - m) - 3E{m)/K{m)]'^' 

as = -Mo + 1.4(a/7)ci[(2-m)-3E(m)/ir(m)]/c2, 
ao = 1.4(r/7)ci[(2-m)-3E(m)/ir(m)]/c2, 



(22) 

and the constants Cq, Ci, C2 were defined in eq. (|18]). 

A steady state providing for equilibrium between the gain and loss is a fixed point 
(FP) of eq. (pip, which is attained when r.h.s. of the equation vanishes. This is possible 
either if s is a root of the cubic equation 

+ a2S^ - s - ao = 0, (23) 

or if s — 1/s — mq = 0, i.e.. 



2 ^u, ± ^ul + 4j . (24) 

The FPs ( p^ lead to a trivial result, as it immediately follows from eq. ( [T3| ) that they 
corresponds to a constant solution with 77 = 0. Therefore, we now focus on eq. (PB)). In 
fact, precisely the same equation ( |2BD can also be obtained from the solvability condition 
of a standard asymptotic perturbation analysis of the full system (we do not present it 
here) . 

Physical roots of eq. (^) are those which are not only real but also provide for rf > Q 
through eq. (0). Noting that the combination (2 — m) — 3E{m)/K{m) in eq. (|19]) is 
negative if m < 0.96, and positive if m > 0.96, we conclude that physical roots must 
belong to the following intervals, in order to guarantee 77^ > 0: 



s < |(mo - Vmo + 4), or < s < |(uo + + 4), if m < 0.96; ^^5) 
|(mo - \/ul + A) < s < 0, or s > |(mo + ^Ju^ + A), if m > 0.96. 

In particular, for uo = 0, physical roots lie in the intervals s < — 1 or < s < 1 if 
m < 0.96, and -l<s<Oors>lifm> 0.96. 

Further, we define q = —1/3 — (aa/S)^ and r = (1/6) (3ao — aa) for eq. (p^). Then, 
as is well known, roots of the cubic equation can be classified as follows: there are one 
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real root and a pair of complex conjugate ones if + > 0, and three real roots if 
g,3 _j_ ^2 ^ g Tii^g^ generally speaking, there might exist up to three physical roots to 
eq. (|23|) . However, a detailed numerical analysis based on the above formulas has only 
revealed cases with two physical roots, or none. In particular, if there is a single real root 
of the cubic equation (|23|), this root is always unphysical, as it does not satisfy the above- 
mentioned condition rf > 0. Two other roots are physical if they are real. Moreover, if 
m < 0.96, both physical roots are negative, s < {uq — \Juq + 4)/2, and if m > 0.96, both of 

them are positive, s > {uq + \Juq + A)/2. Because s is the velocity of the traveling CnWs, 
in these two cases we will call the waves left-moving and right-moving ones, respectively. 

Equation (|2T|) also provides for a necessary, although, strictly speaking, not sufficient, 
information about the stability of the CnW solutions. Indeed, a necessary condition is 
that the corresponding FP be stable within the framework of eq. (|2T|) . It can be checked 
that the sign of r.h.s. of eq. (0) always changes from positive to negative when passing 
(from left to right) through FP corresponding to a larger value of the speed |s| (it can 
be demonstrated that it always corresponds to a larger value of the amplitude r] too), so 
this FP is obviously stable as a solution to eq. (pT]). Similarly, FP corresponding to the 
smaller values of |s| and t] is always unstable, playing the role of a separatrix bordering 
the attraction domain of the stable FP. Direct simulations (see below) have demonstrated 
that, whenever the CnW solution corresponds to FP is stable as the solution to eq. (plf), 
it turns out to be as well stable as a periodic-wave solution of the full system of eqs. (0) 
and d). 

It is relevant to mention that quite a similar conclusion was conjectured analytically 
and proved by means of direct simulations in ref. for SW solutions to eqs. (0) and 
(I). The SW solutions also appeared in pairs, the pulses with larger and smaller values 
of the amplitude being stable and unstable, respectively. 

Note that all the above conditions, for instance, eqs. (^), involve only the mean value 
Uq of the field in the main {u-) subsystem, but not the mean value Vq of the extra field. 
Therefore, we conjecture that the stability of the CnW solutions does not depend on vq. 
This point will be confirmed by numerical simulations reported in the next section. Thus, 
for fixed values of the gain and loss parameters a, 7 and F, the stability of the periodic 
waves is determined by uq and m, the latter constant being related to the period L, see 
eqs. (H). This suggests to display stability regions for the CnW solutions in the (L, uq) 
parametric plane, if the values of a, 7 and F are fixed, which will be done below. 

Adopting the principle stated above, i.e., that the CnW solution with the larger values 
of |s| and 77 is stable, we can summarize predictions for the stability in the form of 
parametric regions; in the next section, the predictions will be checked against direct 
simulations of eqs. (|^) and (^. Figure 1 shows the stability region predicted in the 
plane {L,uo) for the case a = 0.15, 7 = 0.05, and F = 0.2. Three different areas can be 
identified in this case: the right-moving waves are stable in the area I, constant solutions 
are stable in the area II, and the left-moving waves are stable in the area III. It can be 
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seen from Fig. 1 that, within the range considered, the left-moving waves disappear if the 
mean value Uq exceeds a certain critical value 0.6), whereas the right-moving waves 
only exist when uq is Zander than another critical value (~ —0.12). 

Note that, when uq is fixed to be zero. Fig. 1 has gaps on the axis of L inside which 
stable waves are absent. Indeed, the left-moving waves exist if 4.0 < L < 6 or L > 8.0, 
and the right-moving ones exist if L > 16.75, i.e., the gaps are L < 4.0 and 6.0 < L < 8.0. 
This situation can be explained in detail in terms of eq.(^). Indeed, it is easy to check 
that, if Mo = 0, two physical roots to eq. (p3D exist if m < 0.7052 or m > 0.999858, 
the nonexistence gap 6.0 < L < 8.0 in Fig. 1 (for uq = 0) being a consequence of the 
discontinuity in the values of m at which eq. (p3D has two physical roots. More accurately, 
m = 0.7052 corresponds to L = 5.997. It is the maximum value of L for which the left- 
moving periodic waves with a single maximum per period can be found; such waves may 
be naturally called single-hump ones. Left-traveling waves reappear at L > 8.0, this time 
having = 2, 3, ... identical humps inside one period (so that the fundamental period of 
the wave is L/N rather than A^). Single-hump right-moving waves appear when L > 16.75, 
which corresponds to the above-mentioned value of m = 0.999858. At this point, a pair 
of complex conjugate roots of eq. (^) bifurcate into two physical ones. 

On the other hand, eqs. (|23|) and (|2^ ) predict that, as m — 0, the period L takes 
the minimum value, Lmin = 3.965, at which stable CnWs may exist. This explains the 
presence of the other above-mentioned gap, L < 4. 

The model easily gives rise to a multistability, i.e., coexistence of several attractors at 
the same values of parameters, an example being overlapping between the areas I and III 
in Fig. 1. The multistability is well confirmed by numerical simulations reported in the 
next section. 

Figure 2 displays stability areas for another set of parameter values: a = 0.1, 7 = 0.05, 
r = 0.15. In this case, there is no overlapping between the areas I and III, which are 
separated by the area II, where a constant solution is the single attractor. Moreover, 
no stable left- or right-moving waves can be found for mq = in this case, which is a 
drastic difference from the situation displayed in Fig. 1. Left-moving waves exist if uq is 
smaller than —0.15, while the right-moving waves cease to exist if Uq < 0.2 and L < 35. 
In particular, when uq = —0.2, the single-hump stable periodic wave exists in an interval 
5.9 < L < 7.2, and when mq = —0.3, the interval is 5.5 < L < 8.5. For these values 
of Uq, multistability takes place at large values of L, not as the coexistence of left-and 
right-traveling waves, but rather as coexistence of stable waves traveling in one direction 
but having different numbers of humps per period. As it was shown in ref. the model 



also supports a stable SW at the same values of the parameters. Therefore, a conclusion 
is that, as L — > cxD, the stable single-hump right-moving wave carries over into the stable 
SW, with the corresponding mean value Uq vanishing in this limit. 

54. Numerical Simulations 
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In this section, we aim to verify predictions for the existence and stabihty of periodic 
waves by direct simulations. To this end, eqs. @ and were integrated by an imphcit 
Fourier pseudo-spectral method (the same as the one used in ref. p2|), with two different 
types of initial conditions, sm{kx) and sech^(fi;x). 

Firstly, for the same set of parameter values, a = 0.15, 7 = 0.05, F = 0.2, as was used 
in Fig. 1, and Mq = 0, various numerical experiments were conducted with different values 
of L and different initial conditions. As a result, stable periodic left-moving waves with a 
single hump per period were found in an interval 4.0 < L < 6.5, at it was also concluded 
that stable waves are absent in an interval 6.5 < L < 7.9. Left-moving periodic waves 
with two identical humps per period were found for L > 7.9. These findings are in very 
good agreement with the predictions of the perturbation theory which were explained in 
detail in the previous section and summarized in Fig. 1. On the other hand, single- hump 
right-moving periodic waves were observed for L > 18.0, which should be compared to 
the above-mentioned theoretical prediction L > 16.75. 

Figure 3 shows profiles of the stable waves found from the simulations for L = 6.0 and 
L = 10.0, with one and two humps per period, respectively. These profiles were obtained 
from the sin(A;x) initial conditions. At these values of the parameters, all the initial 
profiles either evolve into a single wave, or decay into a trivial state. It can be checked 
that the numerically found values of the amplitude and velocity of the wave agree well 
with the analytical predictions, see caption to Fig. 3. 

It is also necessary to check the prediction that the mean value Vq of the extra field 
does not affect the stability of the periodic waves. To this end, we conducted a numerical 
experiment, fixing L = 12 and uq to be 2.0 or —1.0. These values correspond to points 
in the areas I and III, respectively, in Fig. 1. Then, we investigated the established wave 
profiles for various values of Vq. A typical result is displayed in Fig. 4, clearly showing 
that the mean value Vq affects neither the stability nor the shape of the -u-component. 

Another important issue to be checked is the multistability. Based on the previous 
results, it is clear that the separation between local maxima in a stable periodic left- 
moving wavetrain (for the same values of parameters as in Figs. 3 and 4) is between 4.0 
and 6.0, as predicted analytically, or between 4.0 and 6.5 as found numerically, if Uq = 0. 
Therefore, when the period is L = 18.0, the number N of humps per period must belong 
to the interval 18/6.5 < N < 18/4.0, which yields = 3 or = 4. Thus, the coexistence 
of two left-traveling waves with the fundamental periods L/3 and L/4 is expected in this 
case. As it was mentioned above, the stability of the periodic waves with the fundamental 
period essentially smaller than L implies that waves in the present model may also be 
stable against perturbations breaking their periodicity. In fact, taking larger values of L, 
we were able to observe, at various values of the parameters, stable periodic waves with 
the fundamental period as small as L/16. 

The single-hump right-moving stable wave was predicted to exist too at the same 
values of the parameters for L > 18.0. All these predictions were confirmed by numerical 
simulations. Stable left-moving waves with three and four humps per period, and stable 
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single-hump right-moving wave are shown in Fig. 5 for L = 18.0. Again, the numerical 
results for the amplitudes and velocities of the waves are found to agree well with the 
theoretical predictions. We stress that no other waves but the ones mentioned here had 
been turned up by simulations at these values of the parameters. 

Lastly, we display numerical results for another set of parameters: a = 0.1, 7 = 0.05, 
and r = 0.15. When uq = 0, no left- or right-traveling stable periodic wave can be 
found, various initial conditions relaxing into a trivial uniform state in this case. Recall 
that precisely this feature was predicted in the previous section (see also Fig. 2). With 
uq = —0.2, stable single-hump left-moving waves are found for 5.8 < L < 7.2 (the ana- 
lytical prediction being 5.9 < L < 7.2); when uq = —0.3, stable single-hump left-moving 
waves are found for 5.5 < L < 8.4 (the analytical prediction is 5.5 < L < 8.5). As is seen, 
the numerical results agree with the predictions very well. In other words, the inter-pulse 
separation for a stable periodic left-traveling wavetrain is, in this case, between 5.8 and 
7.2 for Uq = —0.2, and between 5.5 and 8.4 for uq = —0.3. Therefore, if the period L is 
large, the multistability must manifest itself as the coexistence of several different waves 
with multiple humps. Figure 6 shows two- and three-hump left-moving waves for fixed 
values of L = 24.0 and uq = —0.3. In this case, the possible number of the humps per 
period may be 24/8.4 < N < 24/5.5, implying = 3 or iV = 4, which is confirmed 
by Fig. 6. On the other hand, when uq becomes positive, there may exist right-moving 
single-hump waves. Several right-moving waves are shown in Fig. 7 for different values 
of Uq and L. 



5. The Zero-Dispersion Limit 



All the results presented above were obtained for the case when the linear gain and loss 
terms in the underlying equations (0) and (^ matter of fact, small perturbations 

in comparison with the dispersive term (third derivative) in eq. (|^). It is also interesting 
to consider the opposite limit of the zero-dispersion system, when the KS-KdV equation 
(0) is replaced by the KS equation proper, so that the system takes the form 

ut + uu^ -v^ = -au^^ - 'yu^xxx, (26) 
vt-u^ = Tv^x, (27) 



[following the arguments presented above, we set c = in eq. (pT])]. 

First of all, it is relevant to mention that, on the contrary to the model based on the 
KS-KdV (Benney) equation, the zero- dispersion model does not support stable SWs. Our 
numerical simulations demonstrate that, with the decrease of the coefficient D in front of 
the third derivative in the KS-KdV equation [in eq. (^ this coefficient was 1], there is a 
finite critical value Dcr, such that no stable pulses can be found for D < Dcr- However, in 
the coupled system (^) and (pTj), stable periodic waves continue to exist, both with one 
and several humps, and, unlike the KS equation proper (see, refs. |2B|, and references 



therein), the present system is not prone to generate dynamical chaos. 
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Typical examples of stable periodic-wave solutions to the system and ( P7| ) with 
one and two humps are displayed in Fig. 8. Quite similar results were obtained at many 
other values of the parameters, while spatiotemporal chaos has never been observed, al- 
though the search for chaos in the system's parameter space was not exhaustive. 

§6. Conclusion 

In this work, we have studied cnoidal-like waves and their stability in the system 
of linearly coupled KS-KdV equations. We first constructed a family of cnoidal-wave 
solutions for the zero-order system, which does not include the gain and loss terms. Then, 
treating the gain and dissipation as small perturbations, and making use of the balance 
equation for the net field momentum, we have found necessary conditions for the existence 
of stable cnoidal waves. Finally, the predictions were checked against direct numerical 
simulations, showing very good quantitative agreement in all the cases considered. The 
obtained results can be summarized as follows: 

1. For mean value vo of the field in the linear subsystem does not affect stability of the 
periodic waves. The stability is determined by the mean value Uo of the field in the 
main subsystem and the wave's period L. 

2. For fixed values of the model's parameters, three different stability regions are found 
in the (L, uq) plane, viz., region I for right-moving waves, region II for trivial 
solutions and region III for left-moving waves. 

3. The three stability regions can overlap with each other, which leads to multistability. 

4. At large values of L, the multistability manifests itself also as coexistence of stable 
waves moving in one direction but having different fundamental periods L/N, with 
an integer A^. In fact, may take large values, which implies that the periodic 
waves are also stable against perturbations breaking their periodicity. 
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Figure Captions 



Fig. 1. Stability areas in the parametric plane {L,uq) for a = 0.15, 7 = 0.05 and 
r = 0.2, as predicted by the perturbation theory. Right-moving cnoidal waves are stable 
in the region I, which lies above the solid curve, left-moving cnoidal waves are stable in 
the region III, which is below the dashed curve, and the constant solutions are the only 
stable solutions predicted in the region II. 

Fig. 2. The same as in Fig. 2 for a = 0.1, 7 = 0.05 and F = 0.15. 

Fig. 3. Numerically found profiles of stable periodic waves for a = 0.15, 7 = 0.05, 
F = 0.2 and Uq = (one period is displayed). The solid and dashed curves show the 
u- and f-components: (a) L = 6.0; (b) L = 10.0. In these two cases, numerically found 
amplitudes of the u - component are, respectively, Au = 3.46 and = 5.85, while the 
perturbation theory predicts A^ = 3.65 and A^ = 5.95. 

Fig. 4. Numerically found stable waves for a = 0.15, 7 = 0.05, F = 0.2, L = 12.0 
and different values of vq. The solid curve is the profile of the u-component, while the 
dashed and dashed-dotted lines are profiles of the f-component obtained for different Vq 
(the M-component is shown by a single curve, as the change of vq produces no visible 
variation in it): (a) uq = 2.0, with f = (dashed line) or vq = 2.0 (dot-dashed line); (b) 
Mo = —1.0, with Vq = —1.0 (dashed line) or t>o = (dot-dashed line). 

Fig. 5. Multistability of stable periodic waves for a = 0.15, 7 = 0.05, F = 0.2, 
Mo = 0, and L = 18.0. The solid and dashed curves represent the profiles of the u- and 
w-components. Panels (a) and (b) show the left-moving waves with 3 and 4 maxima per 
period, and panel (c) shows the right-moving wave with the single maximum. 

Fig. 6. The same as in Figs. 5(a) and 5(b) for a = 0.1, 7 = 0.05, F = 0.15 and 
L = 24.0. 

Fig. 7. Stable right-moving cnoidal waves for a = 0.1, 7 = 0.05 and F = 0.15. The 
solid and dashed curves show the u- and v - components for (a) L = 12.0, uo = 2.0, (b) 
L = 15.0, Mo = 1.0, (c) L = 30.0, Uq = 0.4. The numerically found and analytically pre- 
dicted values of the amplitude of the u - component in these three cases are, respectively, 
A^ = 5.48 and 5.70; A^ = 4.50 and 4.68; A^ = 5.23 and 5.42. 

Fig. 8. Typical examples of stable periodic waves with one (a) and two (b) humps per 
period, found by means of direct simulations of the zero-dispersion model based on Eqs. 
(H) and (P), for a = 0.15, 7 = 0.05, F = 0.2, and L = 10 (a) or L = 11 (b). In both 
cases, the mean value of the w-field is Uq = —2. The profiles of the u - field are shown. 
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